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Abstract. All observations by the aperture photometer (PHT-P) and the far-infrared (FIR) camera section (PHT- 
C) of ISOPHOT included reference measurements against stable internal fine calibration sources (FCS) to correct 
for temporal drifts in detector responsivities. The FCSs were absolutely calibrated in-orbit against stars, asteroids 
and planets, covering wavelengths from 3.2 to 240 fim. We present the calibration concept for point sources 
within a flux-range from 60 mjy up to 4500 Jy for staring and raster observations in standard configurations and 
discuss the requisite measurements and the uncertainties involved. In this process we correct for instrumental 
effects like nonlinearities, signal transients, time variable dark current, misalignments and diffraction effects. A 
set of formulae is developed that describes the calibration from signal level to flux densities. The scatter of 10 
to 20 % of the individual data points around the derived calibration relations is a measure of the consistency 
and typical accuracy of the calibration. The reproducibility over longer periods of time is better than 10%. The 
calibration tables and algorithms have been implemented in the final versions of the software for oflline processing 
and interactive analysis. 
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1. Introduction 

In its broadest meaning, we define calibration as trans- 
forming the specific, nowadays digital, output of a sci- 
entific instrument to physical units. In our case these are 
flux densities at given wavelengths and sky positions. This 
transformation generally varies with every change of the 
instrument set-up, hence the complexity of the calibration 
task increases with the number of instrument configura- 
tions used. Initially it is derived on the grounds of the 
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known instrument geometry and the relevant optical and 
electrical properties, that as a whole we will refer to as 
the ideal instrument model. Subsequently this model is 
refined and becomes more empirical, in order to match 
the measured data. 

The task of calibration can be divided into three parts: 
First, development of the instrument model and deter- 
mination of the instrument parameters that are assumed 
to be unchanging and that can be measured in the lab- 
oratory, e.g. filter transmissions, aperture diameters, etc. 
Second, the determination of the open parameters of the 
ideal instrument model that can be determined only in- 
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situ, i.e. with the instrument built into the satellite and 
in the real space-environment. And third, the determina- 
tion of deviations from the ideal instrument model that 
cannot be removed by adjusting parameters, and require 
a new functional property. These modifications to the in- 
strument model are found empirically, and generally origi- 
nate in simplifications. Because open parameters and non- 
ideal instrumental effects are usually intertwined, an iter- 
ative process is required to separate and quantify all the 
contributing effects. 

This paper presents the photometric calibration of 
staring mode- and simple raster mode observatiouf^ of 
point sources with the P- and C-sections of ISOPHOT 
(Lemke et al. 1996| ), which is one of four scientific instru- 
ments on board ESA's Infrared (IR) Space Observatory 
ISO (Kessler et al. |l99|). Unlike CCD devices, the de- 
tectors for the Mid-IR (MIR) and Far-IR (FIR) are far 
less stable and exhibit a continuously changing relation 
between signal and incident flux. Thus the two Fine 
Cahbration Sources (FCS) built into ISOPHOT played 
a crucial role as stable references for the photometry and 
most of this paper will describe their empirical calibration. 

The data collected for this task represent the largest 
part of all specific calibration observations during the mis- 
sion. This resulted in a fairly homogeneous block of data. 
Its analysis and comparison to modeled spectral energy 
distributions (SEDs), of the observed celestial standard 
sources yielded most of the results presented here and 
drove several refinements to the ideal instrument model. 
An additional difficulty was the large number of possible 
instrument configurations, which was limited somewhat 
by considering only one standard aperture for each filter 
band of the aperture photometer (see Table ||). The set of 
instrument configurations and modes we treat herein de- 
fines a well-understood baseline within the large parame- 
ter space, where absolute calibration errors are expected 
to be minimal. The calibration of further configurations 
and modes, like chopped observations, extended source 
photometry, or non-standard apertures, is left to future 
pubHcations. 

We start with a brief review of the instrumental de- 
sign with some emphasis on the internal reference sources 
(Sect. H), followed by an outline of the calibration strat- 
egy (Sect. 1^). Sect. (^ continues with a description of 
the corrections applicable to the detector signal. Sect. (||) 
presents the celestial calibrators and Sect. (^) describes 
those corrections imposed by photometric constraints. We 
derive the FCS calibration tables in Sect. (^ and present 
the final flux calibration, with its mathematical descrip- 
tion and a discussion on accuracy and reproducibility, in 
Sect. (||). A summary constitutes Sect. (||). 



2. Instrument Design 



O Sky 




Fig. 



Detector 



1. Instrument schematic from a calibrator's point of 
view. The detector compares IR radiation directly from the 
sky and from the internal reference source. Stability of the 
detector is required only for the time interval of the two 
measurements. 



2.1. Optical System 

A schematic instrumental set-up is shown in Fig. |^. A 
more detailed description can be found in Klaas et al. 
( 1994 ) or Laureijs et al. ( ^001 ). The set-up allows a choice 
between two sources of IR radiation via a chopper mirror: 
i) the telescope and ii) the internal FCS. The radiation 
of either source is analysed by an optical system through 
filters, apertures, etc., and reaches finally the IR detectors. 

The telescope is a Ritchey-Chretien design, diffraction- 
limited at 5 /im, with an entrance aperture 60 cm in 
diameter and an /-ratio of 15 (Kessler et al 



1996) 



Excluding the AOT P32. 



Radiation reaches the chopper mirror within ISOPHOT 
via a pyramidal mirror, which is centred between the 4 
ISO-instruments. The sky is projected onto the focal plane 
with a scale of 0.04363 mm/". 

2.2. Detectors 

The subsystems PHT-P and PHT-C are a multi-aperture 
photometer with 3 single pixel detectors and a FIR camera 
comprising 2 detector arrays of 3 x 3 and 2x2 pixels. 
The detectors of PHT-P are PI (3.15-17.5 ^im), P2 (15- 
30 /im), and P3 (40-130 ^m) with 10, 2, and 2 filter bands 
respectively. Each detector can be combined freely with 
one of 11 circular field apertures with diameters between 
5"and 180". Only one filter band is measured at a time. 
The PHT-C pixel sizes are 1.9x1.9 mm and 3.9x3.9 mm 
with gaps of 0.1 mm between pixels. The projected pixel 
sizes on the sky are 43.5" x 43.5" for ClOO and 89.4" x 89.4" 
for C200. 

The detectors, made from doped silicon or doped ger- 
manium, are extrinsic photoconductors operated at tem- 
peratures of 1.7 to 3.6 K, depending on material. Their re- 
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sponse to IR radiation is characterised by the detector re- 
sponsivity, i?, which is the parameter most relevant to cal- 
ibration. It is given by the relation between detector-signal 

and incident in-band power R = -^^ — where S is 

the detector signal expressed in V/s for an integrating am- 
plifier. S'dark is the relaxed signal under dark conditions, 
Cint is the integrating capacity (140 fF for C200, 90 fF for 
others) and P represents the in-band power in Watts. 

The responsivity of these detectors is not constant. 
Depending on material, variations of a factor of 3 during 
one orbit were observed. Changes of R are induced by the 
ionizing radiation in space and by the previous sequence of 
IR fluxes to which the detector has been exposed. We con- 
sider 7? to be stable over typical time intervals of around 10 
to 20 min. To bring R back to a nominal value, the detec- 
tors routinely underwent a curing-procedure after having 
crossed the Radiation Belts during perigee and before reg- 
ular observations started ("science window") (see Laureijs 
et al. 1996 2001). These procedures consisted of detector 
specific combinations of detector-heating, bias increase or 
illuminator fiashes (Lemke at al. 1996) and brought the 
responsivity back to within 5 % of its nominal value. The 
doped-Ge detectors, which showed the biggest responsiv- 
ity changes, underwent a second curing procedure after 
8 hours, close to apogee. 

2.3. Fine Calibration Source (FCS) 

The instrument hosts two FCSs, of which FCS 1 was 
used for the staring and raster observations and, hence, 
is chiefly treated in this paper. Both devices were man- 
ufactured at MPIA Heidelberg. Each FCS contains two 
thermal radiation sources (TRS) made of small diamond- 
plates of 1 X 1 X 0.1 mm, coated with Ni-Cr. The plates 
are heated electrically with stabilised heating powers up 
to 50 mW, providing highly reproducible IR fluxes. The 
stabilisation circuit works digitally, with a resolution of 
12.21 yuW. This assures a flux stability of better than 4% 
for heating powers above h — 120 fj,W, and 2% above 
h = 264 /iW. 

The radiation of both sources is combined in a dia- 
mond beam-splitter. The beam of TRS 1 is reflected with 
practically no losses, while the beam of TRS 2 is trans- 
mitted through the material and attenuated by a factor of 
~1000. Only one TRS is heated at a time. The attenuated 
TRS 2 is operated with the doped-Ge detectors P3, ClOO 
and C200, which show higher responsivities, whereas the 
doped-Si detectors PI and P2 are used with TRS 1 in the 
stronger reflected beam. Mirrors focus the emerging ra- 
diation in such a way that a beam with the same aspect 
ratio as the telescope is emitted, while the TRSs are im- 
aged onto the chopper mirror, making the FCS appear to 
the detectors as a homogeneously illuminated source. 

The in-band power provided by the FCS could be 
adjusted by changing its heating power. The fundamen- 
tal calibration curves, in-band power versus FCS heating 
power, were initially determined by a grey body model 



(Schulz 1993), fitted to a few already known data points. 
The grey body is defined as eB^{T), which is the Planck 
function multiplied by a wavelength-independent emissiv- 
ity e. The temperature, T, is linked to the FCS heating 
power, ft,, by an empirical function, T{h) = a * -I- 7 -I- 
(5atan(P/2), where 7 is the temperature of the optical 
support structure (2.76 K), and a, (3 and 5 are dimen- 
sionless constants. Typical values for these constants of 
FCS1/TRS2 arc 25.8 ± 1.8, 0.518 ±0.006 and 0.07 ±0.04 
respectively. An attenuation factor which modifies the ef- 
fective solid angle individually for each detector subsys- 
tem was introduced, describing flux losses that were not 
predicted by the grey body model. 

The advantage of the grey body model was to allow 
predictions for other filter bands, when only a few data 
points per detector were available. To increase accuracy, 
once enough data points were available, the covered ranges 
were interpolated by smooth low-order polynomials and 
the model was used for extrapolations only. 

3. Calibration Strategy 

Given the drifting detector responsivity, the calibration 
had to be reestablished periodically by observing a known 
reference source. To avoid time-consuming slews across the 
sky to celestial calibration sources, we used the FCSs as 
an intermediate reference that could easily be compared 
with the telescope beam (Fig. |l|). AU ISOPHOT observa- 
tions that were implemented as Astron omical Observing 
Templates (AOT)(see Klaas et al. |l994| ), were designed to 
contain at least one FCS measurement. 

The absolute calibration of these secondary standards 
was established in-orbit by comparing their output to 
known celestial standards, an activity that started during 
the performance verification phase (PV), but continued 
throughout the mission due to the visibility constraints of 
some of the sources. To minimize errors due to potential 
detector nonlinearities, the FCS heating power was ad- 
justed so that the emitted FCS flux roughly matched the 
flux emitted by celestial source and background. The lim- 
ited linearity of the system demanded celestial calibration 
sources at all flux levels in all 25 filter bands. 



4. Deriving the Detector Signal 

To calculate the detector responsivity, i?, we need to deter- 
mine the detector signal. Ideally the signal is proportional 
to the detector current, exhibits only Gaussian distributed 
noise, and remains constant over periods of constant de- 
tector illumination. To approach this ideal in practice, a 
number of effects must be treated, which lead to the signal 
conditioning procedures described below. 



4.1. Integration Ramps 

Integrating cold (« 3 
(Dierickx et al. 
rents of typically 10 



1989|) 



-16 



K) readout electronics (CRE) 
are used to amplify the cur- 
. 10" ^■^ A that flow through the 
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Fig. 2. Ramp nonlinearity for pixel 8 of CI 00. Deviations 
of many integration ramps from a straight line are plot- 
ted against CRE output voltage. This relation is used for 
ramp-linearisation. 



ISOPHOT photoconductors (Lemke et al. |1996|). The pho- 
tocurrent is measured from the rate at which the voltage 
at the charge capacitor within the CRE increases with 
time. The voltage is sampled at regular time intervals for 
a given duration before it is reset (non-destructive read- 
outs, NDR, and destructive readouts, DR). The readouts 
of all channels of a CRE unit are time-multiplexed and 
sent via a single line to the external electronics unit (EEU) 
outside the cryostat. There they are further amplified and 
digitized by a 12-bit analog-to-digital-converter (ADC). 
The sequence of samples between two resets is called an 
integration ramp and ideally fits a straight line. In the 
following we will refer to the slope of the fitted line as 
"signal", measured in V/s. Multiplication by the integra- 
tion capacity, for which the design value of 90 fF (for C200 
140 fF) is adopted, leads to the photocurrent in Amps. 

4.2. Nonlinearities of Integration Ramps 

The actual integration ramps arc, however, not perfectly 
straight (see Fig. |) (Schulz [19931) . Two effects can be sep- 
arated: i) The AC-coupled CRE circuit is not a perfect 
integrator. The output voltage can rise only to the level of 
the detector bias times the small amplification of the cir- 
cuit (wlO). This saturation level is reached asymptotically 
via a typical RC-loading curve. High bias voltages above 
10 V have sufficiently high saturation levels («100 V and 
more) that they show practically linear integration ramps 
within the CRE output range of «2 V. However, for the 
smaller biases of P3, ClOO and C200, the saturation volt- 
age drops such that the curvature of the loading curve ap- 
pears within the dynamic range of the integration ramp 
(debiasing). ii) Secondly, all integration ramps show devi- 
ations from a straight line that always appear at the same 
CRE output voltage, regardless of the level from which 
the ramp started to integrate. Due to the direct link to 
the CRE output voltage we attribute this component to 
an intrinsic nonlinearity of the amplifier within the circuit. 



The correction algorithm is based on the observation 
that all integration ramps of one detector pixel can be 
matched by conserving the CRE output voltage of the in- 
dividual readouts, while linearly stretching and shifting 
the time- axis of each ramp (Schulz [1993 ). In this way an 
average integration ramp was determined for each detec- 
tor pixel from a number of measurements with ramps of 
large dynamic range. The difference between a straight- 
line fit to this ramp and the ramp itself (Fig. ^ is used to 
linearise individual ramps. The maximum absolute values 
range between 40 and 100 mV over the maximum CRE 
dynamic range of 2 V. The correction depends only on 
CRE output voltage, and has no further dependences on 
readout frequency, in-band power or time. 

The reliability of the correction is deduced from the 
standard deviation of the individual ramps from the av- 
erage integration ramp. Over most of the dynamic range 
it is below ±10 mV, with a tendency towards larger scat- 
ter for longer wavelength detectors. The highest standard 
deviation appears in the central pixel of ClOO, showing 
±20 mV. The higher scatter towards long wavelength de- 
tectors with smaller biases is interpreted as resulting from 
a simplification in the correction algorithm. Relating the 
corrections only to the CRE output voltage neglects the 
fact that the curvature due to debiasing also depends on 
the reset level of the integration ramp, which shifts slightly 
w.r.t. the CRE output voltage. Hence an additional scat- 
ter appears, increasing with smaller bias. Considering the 
maximum standard deviation of ±20 mV and an aver- 
age dynamic range of 1 V (half the maximum dynamic 
range), we estimate the typical residual systematic error 
to be about ±2 % after correction. 



4.3. Ionizing Radiation Effects and Deglitching 

The output signal is disturbed by energetic particle hits 
(glitches), mostly protons and electrons. A typical distri- 
bution of fitted slopes after ramp linearisation is shown 
in Fig. ^ The asymmetry of the distribution is caused by 
glitches that instantaneously increase the charge on the ca- 
pacitor during the integration, so that the average slope 
of the affected integration ramp is increased. Typically, 
the readouts following the glitch continue integrating as 
before; however, stronger hits can affect the detector, lead- 
ing to long-lasting responsivity changes. A number of al- 
gorithms have been developed to remove signals affected 
by ghtches (Gabriel et al. |1997D. 



4.4. Digitization Effects 

All deglitching algorithms still leave an asymmetry in the 
signal distribution. For stronger asymmetries the simple 
average is often quite far from the peak of the distribu- 
tion. The median comes closer to the peak; however, it 
is a bad choice for very weak signals and short integra- 
tion ramps, where the bins of the A/D-conversion become 
significant. The AC-converter covers the range from —10 
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Fig. 3. Histogram showing the signal distribution in a long 
measurement of the central pixel of ClOO. The tail towards 
higher signals is caused by high energy particle hits and can 
be removed by appropriate algorithms. 



to +9.995 V. Dividing by the gain factor of the analog 
electronics and 2^^ for 12 bits, the least significant bit is 
equivalent to 6.1 mV. For integration ramps of 1 second 
this would yield a detector current of 5.5 x 10^^^ A. As 
an example, the 100-/im channel would produce discrete 
slope levels in about 20 mJy intervals. Since the result 
of the median is always an existing value of the sample, 
the accuracy is limited to the separation of the intervals. 
In this case, better results can be obtained by fitting a 
Gaussian to the signal distribution. 



The simplest solution is to determine the level when the 
signal has practically reached its final value and no trend 
is detected within the noise. To filter out the cases that 
need special treatment, a 'Stability Check', implemented 
in the PHT Interactive Analysis (PIA) software (Gabriel 
et al. 1997), is conducted. The method tests successively 



smaller pieces of data towards the end of a measurement 
for stability at a 95 % confidence level. The smallest piece 
is limited to either 7 data points or 8 sec. Stable signal lev- 
els were detected for PI, P2 and ClOO for about 65% of 
our calibration measurements. Typical measurement times 
were between 64 and 256 sec. For fainter signals, however, 
the relaxation times generally increase and other methods 
had to be applied to predict the final signal level from only 
the first part of the transient. 

4.5.2. Empirical Transient Fit 

Transients converge asymptotically towards a final value, 
Soo- Empirical functions that fit this characteristic are 
usually constructed from exponentials. Since long term 
detector drifts bound the integration time above, typical 
durations of 'staring' measurements at one flux level range 
between 64 and 256 sec. We use an empirical function of 
the form S{t) = S'oo + CS'o — 5'oo)e~'^*. 5*0 is the signal before 
the flux change, t is time and c is a free parameter. Since 
in most cases Sq is not known, it is also left as a free pa- 
rameter to fit. Stronger transients that are 'remembered' 
by the detector during successive measurements can be 
eliminated by fltting a 'baseline' of the same expression, 
if more than one measurement of the same flux level is 
obtained (see Fig. 0). 



4.5. Detector Transients and Long-Term Drifts 

After a change of IR flux, the photocurrent shows a 
characteristic behaviour before it approaches the flnal 
value asymptotically, which we call hereafter 'transient'. 
Relaxation times of several minutes are common. The sec- 
ond effect, which we refer to as drift, is also observed as 
a gradual variation of the responsivity; however, it is on 
timescales of hours rather than minutes. Drift is due to 
the continuous bombardment of the detector material by 
high energy particles in space (e.g. see Blum et al. 1990] ), 
but can also be triggered by strong flux changes. In the 
following we describe in more detail how the transient cor- 
rection was applied to the FCS calibration data, since this 
was an important step in eventually achieving consistent 
photometric results. 

4.5.1. Solutions for Transients 

The observed detector current is the sum of dark current 
and a current propo rtiona l to th e number of incident IR 
photons (Bratt 1977; Sclar 1984 ). However, this is correct 
only for a system that has reached equilibrium after each 
flux change. To derive signals that are independent of mea- 
surement time, we need to determine the signal as f — > oo. 



o 

> 




bg 
FCS 



initial drifting signal 
corrected signal 
background signal 
internal calibration signal 
source signal 



measurement time [s] 

Fig. 4. Observations in 'nodding mode' with a linear 3- 
step raster used for faint sources. Several changes of the 
pointing between a strong background and a faint source on 
top of it make the source signal appear as a modulation on 
top of a strong transient. This is corrected by dividing by a 
baseline fitted to all background levels. The sequence starts 
and ends with FCS measurements. 
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4.5.3. Automatic Validation of Transient Fits 

Automatic procedures were developed to reduce the cali- 
bration data homogeneously. The routines include stabil- 
ity test, empiric al tra nsient fit and error discrimination^ 
(Huth & Schulz 199^ ). If the stability test cannot locate 
a stable part of the measurement within the criteria given 
above, a transient fit is initiated. The criteria that were 
developed during intensive tests to eliminate dubious re- 
sults from transient fitting are as follows: 

(i) We calculate the average of the last 30 % of the 
measurement and call it Sso- The valid fit result must 
deviate by no more than ±50 % from S30. 

(ii) A straight line is fitted to the last 40 % of the mea- 
surement and extrapolated to a point at 2 times the mea- 
surement time. 80 % of the difference between this point 
and 5*30 is the maximum difference allowed between the 
signal level predicted by the transient fit and 5*30. 

(iii) For each fitted transient, the term (1 — 
Si/ S2)/ {ti — ^2) is evaluated twice, where S and t are the 
slopes and times of either the first two, or the last two, 
data points in the measurement. Absolute differences be- 
tween these two results that are smaller than 0.001 s^"^ 
are rejected. This and criterion (ii) buffer the curvature of 
the fit against small values since small curvatures result 
in very uncertain asymptotes. 

In case these criteria fail, we take 6*30 as the best guess 
for Soo but the result is flagged as being less reliable. This 
applied to about half the cases where transient fits were 
attempted. We caution that the criteria were tuned empir- 
ically to this large, but special, sample of FCS calibration 
measurements. We assume, therefore, that they work best 
for measurement times > 64 s. 



4.5.4. Signal Uncertainties 

In the automatic processing we distinguish three types 
of measurement. The first consists of measurements that 
reach stability within their integration time. Here the sta- 
tistical error is the standard deviation of the mean of the 
stable part of the measurement. The second type consists 
of transient fitted m easurements that meet the c riteria 
above, and we use ^Jl/{n - 1) YJliSf""^^ - 5','°°^'') as er- 
ror, where n is the number of valid slopes in the measure- 
ment, and 5'™°'!' and S™^'^ are the modeled and measured 
slopes respectively. If the transient function were a con- 
stant, i.e. a horizontal line, the expression would be the 
standard deviation. We find that using the above equation, 
instead of the standard deviation of the mean, provides 
more realistic errors. Finally, in cases where no stability is 
found and the transient fit fails, the error assigned is the 
standard deviation of the last 30% of the measurement. 



Note that the quality check function was specifically devel- 
oped to process FCS calibration measurements automatically. 
Although this check is not available in the PIA software, the 
described transient correction can be interactively applied by 
the user, so that our conclusions are not necessarily limited to 
our specific set of measurements only. 
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Fig. 5. The diagram shows the relations of signals of the 
same in-band power measured with different readout tim- 
ing settings. The signals on the y-axis are measured using 
a 1/4- s reset interval and a data reduction factor of 1. 
The reset interval and data reduction factor, respectively, 
for the signals of the x-axis are given by the two numbers 
in each diagram. The data points follow a linear relation, 
with slope different from unity and non-zero x-intercept. 
This example was measured for detector ClOO. 



Thus the errors of the three types also weight the results 
with respect to their systematic uncertainties. 

4.6. Signal Dependence on Readout Timing 

The readout timing (ROT) of the ORE was selected indi- 
vidually for each measurement according to the observer's 
flux estimate, to give a maximum dynamic range without 
saturation and thus a minimum of readout and sampling- 
noise. It is controlled by 3 parameters: the time interval 
between two destructive readouts, trcsot! the number of 
non-destructive readouts, tindr; and the data reduction 
factor. For the range of 128 s to 1/64 s for Preset, depend- 
ing on the detector, a total of 14 to 15 combinations was 
actually used during operations. 

We observed several non-variable^ celestial sources 
(stars, planets) with all ROT combinations suitable for in- 
band powers that cover a range of a factor of 16 above and 
below that actually observed. Care was taken to ensure the 
stability of the detector response during the measurement 
series by repeating the set of ROT combinations in reverse 
order. Ideally one expects to flnd the same signal within 
the errors for the same source, regardless of CRE set-up. 
However, even after correcting for ramp nonlinearity, we 
derived different signals for the same source. The signals 
show a strong dependence on the reset interval, especially 
for those cases where part of the integration ramp exceeds 
the dynamic range and is discarded due to saturation. 



For the time of the observation. 
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To achieve consistency, we correct all signals so they 
appear as if measured with a common reset interval of 
irosct = 1/4 s and a data reduction factor of 1. The sig- 
nal, S, is corrected according to S" = ^0 + ^1 x S, where 
the parameters AO and Al depend on the ROT set-up. 
They were determined from a least-squares fit to the cal- 
ibration data. An example is shown in Fig. ^ It should 
be noted that AO is generally different from zero and the 
slope Al is significantly different from unity. All pixels of 
equal detector-subsystems were found to follow the same 
relation, suggesting that a pure CRE-effect is observed. 
Again the residual errors of the correction are systematic 
in nature. The scatter of the calibration data around the 
fits suggests typical systematic errors of no more than 5 % 
and probably better than 2%. 

4.7. Dark Current 

Even under completely dark conditions and after relax- 
ation of transients due to previous illumination, a detec- 
tor current is measured, arising from thermal generation 
of free charge carriers and surface leakage conductance 
(Bratt 1977). This value must be subtracted to obtain the 
pure photocurrent. The dark current had already been 
determined for each detector pixel in the laboratory; how- 
ever, in orbit higher dark currents than on the ground were 
generally observed. A comparison is given by Lemke et al. 
( |l996| ) . We interpret part of the increase as due to the faint 
end of the glitch distribution, where many small glitches 
produce a continuous current (see Fig. ||). Dark signals 
measured randomly over the mission versus relative or- 
bital position are plotted in Figs ^ and ^ The range to 
1 on the X-axis corresponds to the full 24-hour revolution, 
starting at perigee. The data show a dependence of the 
dark signal on orbital position. The systematic rise along 
the orbit is attributed to high energy particles in space. In 
particular, the steep rise at the end indicates the begin- 
ning of entry into the Radiation Belts. The uncertainty in 
the dark signals is due to detector noise and limitations in 
correcting detector transients. A low order polynomial fit 
to these diagrams is used to model the time-dependence 
of the dark current over the orbit. 



the full observable flux range between 3 and 240 /im. The 
upper limit of the absolute flux range is determined by the 
dynamic range of the CRE and the detector responsivity. 
The lowest signals are mainly determined by the celestial 
background and the detector responsivity. 




A. [[im] 

Fig. 8. The coverage of the flux domain accessible to 
ISOPHOT by the three kinds of celestial standard. 

In Fig. H we illustrate the coverage of the flux domain 
accessible to ISOPHOT by celestial standard sources. The 
largest background components are zodiacal thermal re- 
emission, peaking around 20 /zm, and Galactic cirrus, 
contributing most at the longest wavelengths. Only at 
the shortest wavelengths does the dark current dictate 
the lowest signals. The dotted line shows typical surface 
brightnesses for celestial backgrounds, converted to an 
equivalent flux density of a point source, observed in the 
standard aperture without background. The upper and 
lower absolute flux limits are plotted as dashed lines. Note 
that the limiting sky background for the overlapping wave- 
lengths of PHT-P and PHT-C is different due to their dif- 
ferent aperture areas. Examples of each of the three types 
of calibration source (planets, asteroids, stars) are plotted 
as solid lines. 



5. Celestial Calibrators 

The absolute photometric calibration of ISOPHOT is 
based entirely on point sources, which constitute the best 
known IR flux densities in the sky. The sources were ob- 
served either for PHT-P at the centre of a standard cir- 
cular aperture that is deflned per filter^ (see Table ||) or, 
for the PHT-C array detectors, at the centre of each in- 
dividual pixel. As the detailed radiation spectrum of the 
FCS was unknown and linearity of the detectors could not 
a priori be postulated, the empirical FCS calibration re- 
quired the availability of calibration standards covering 



The standard apertures were chosen to include at least 
90 % of the energy of the point-spread- function (PSF). 



5.1. Stars 

A large database of SEDs of standard stars was prepared 
in support of the ISO calibration. Important aspects of 
the source selection were non- variability, no known stellar 
companions, no IR excess, low sky background and good 
visibility by ISO throughout the mission. 

The ISO Ground Based Preparatory Programme 
(GBPP) (Jourdain & Habing |1992| ; van der Bhek et al. 
[1992 ; Hammersley et al. 1998) provided SEDs between 1 
and 160 fim b y inte rpolating in the Kurucz stellar model 
grids (Kurucz 1993 ), using temperature, surface gravity 
and metallicity. The temperatures were derived either 
from the IR Flux Method (Blackwell et al. |l99ll [l998h 



or the V-K versus Tcff relationship (Di Benedetto 



1993 
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Fig. 6. The evolution of the dark signal with relative orbital position. The data points were taken randomly during the 
mission and plotted versus orbital position, where the range to 1 corresponds to the full 2^ hour orbit, starting at 
perigee. 
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Fig. 7. As for Fig. [| but for the detectors CI 00 and C200. 



199S| ) . Surface gr avitie s and metallicities were from Cayrel 
de Strobel et al. ( 1992 ). The spectral shapes derived from 



the grids were then normalized using near-IR (NIR) pho- 
tometry (Kn or K-band). The attributed errors range be- 
tween 3 % and 5 % which is dominated by the error in con- 
version to an absolute flux density scale. The random scat- 
ter between the measured and predicted fluxes at 1 /im is 
shown to be better than 1.5% (Hammersley et al. 1998 ) 



Another calibration programme by Cohen et al. (1996) 
provided empirical SEDs in the range from 1.2 to 35 /im 
by splicing together measured spectral fragments of cool 
K- and M-giants and calibrating them by absolutely cali- 
brated Kurucz models of Sirius and Vega. The measured 
"composite" SEDs are used to deriv e so-c alled "template" 
SEDs for fainter stars (Cohen et al. 199E ), assuming that 



the intrinsic spectral shape depends only on spectral type 
and luminosity class. The absolute flux level is set by well- 



characterised NIR and MIR photometry, including IRAS 
data. To extend those composite-SEDs not observed as far 
as 35 /im by the Kuiper Airborne Ob servatory, E ngelke 
functions (Engelke |l992| ; Cohen et al. |l995| , | l996|) , with 
effective temperatures from Blackwell et al. ( 1991 ), were 
used. To create FIR extensions to support the longest 
wavelengths of ISOPHOT (240 /im), continuum model at- 
mospheric spectra were attached to the empirical compos- 
ites and extended to 300 /im (Cohen at al. 1996 ). The typ- 
ical absolute accuracy (i.e., in Fa or F^ at any wavelength) 
of the templates is about 3 %, while the FIR extensions 
ranging from 25 to 300 /im have computed uncertainties 
of about 6 %. 

Note that both programmes deliberately use the 
same Vega model SED for their zero point definitions. 
Furthermore, the measured zero points of the various pho- 
tometric systems were determined to ensure that there 
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Table 1. The list of stars, planets and asteroids used for ISOPHOT calibration. The flux densities F^, are given at the 
reference wavelengths of the filters P-11.5 and C-60, depending if the respective standard was used in the range 3-30 fim 
or 45-240 fim. In addition we list the magnitudes in the ISOPHOT photometric system with the Vega model SED as 
zero point. The identifiers MC for M. Cohen and PH for P. Hammersley indicate the SED version that was actually 
used for calibration. Asterisks indicate absolute validation by MSX (Cohen et al. 2001). For solar system objects we 
give the brightness range observable over the ISO mission. 



would be no systematic differe nces depending upon which 
data were used (Cohen et al. 1999 ). Specifically for ISO, 



both programmes supplied electronic versions of their 
SEDs with spectral resolutions of typically 50 to 300. The 
objects are listed in Table 0. 

The irradiances of a number of primary (Sirius, 
Vega), secondary (bright K- and M-giants) and tertiary 
(template-SEDs) calibration stars have been absolutely 
validated by a dedicated radiometric calibration experi- 
ment carried out on the US Midcourse Space Experiment 



(MSX: MiU et al. |l99|) by Cohen et al. (|20M). The 1.2- 
35 jMo. spectra of seven of the stars in Table ^ have been 
validated by this means. We note that because of its ob- 
served deviation from the model SED longward of about 
17 /J,m, we do not use observations of Vega to calibrate 
any long wavelength filter. 

Brighter flux levels in the FIR (i.e. filters with refer- 
ence wavelengths >50 ^m)) had to be calibrated by other 
objects, since even the brightest standard star measured 
by ISOPHOT (a Boo) drops to less than 2 Jy at 200 ^m 
(see Fig. I). 



5.2. Planets 

Uranus and Neptune are commonly used as submillimeter 



and FIR standards (Hildebrand et al. 1985; Orton et al. 
19861; GrifRn & Orton |1993| ). They provided the highest 



calibrated flux levels that ISOPHOT could measure at 
long wavelengths. Mars was already too bright when it 
became visible to ISO. 

In an effort to homogenize the calibration of ISOPHOT 
and the Long Wavelength Spectrometer (LWS) in the long 
wavelength range, for the final version of the offline pro- 
cessing software (OLP 10) we changed the planet models 
used in earlier versions (Schulz et al. |1999D . The ones sup- 
plied by Abbas (1997, priv. comm.) were replaced by mod- 
els produced by GrifRn & Orton for LWS. The models are 
calibra ted w ith 0.35 to 3.333 mm JCMT data (Griffin & 
Orton 1993), which in turn were calibrated against Mars 
(Wright |1976|; Wright & Odenwald |1980[). At the short 



wavelength end, the model was constrained b y Voy ager 
IRIS data from 25 to 50 /imby Hanel et al. (|l986|) for 
Uranus, and Conrath et al. (|198S|) for Neptune. The tem- 
perature structure and composition of H2, He and CH4 
was taken from Voyager radio occultation experiments by 
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Lindal et al. ( |1987| ) for Uranus, and by Lindal ( |1992| ) for 
Neptune. 

We compared the models of Griffin & Orton (GO) with 
the ones of Abbas (AB). Above 45 fim the GO model 
of Neptune is wlO% brighter than AB. In contrast, for 
Uranus the GO model fluxes are about 3 to 10 % lower 
than AB above 50 fim and deviate no more than 15 % 
above 20 /im. At shorter wavelengths the discrepancies 
between AB and GO rise dramatically for both planets. 
However, these are no longer relevant to the PHT cali- 
bration, since the C_50 filter passband starts only around 
40 /im. We take these differences as indicative of how well 
the fluxes from Uranus and Neptune are known at FIR 
wavelengths and conclude that, for our purposes, they are 
still within a ±10% margin. 



5.3. Asteroids 

A small group of asteroids was chosen to be used as FIR 
calibration standards. They populate the intermediate- 
flux level gap that appears at wavelengths from 45 to 
200 fim between Uranus, Neptune and the brightest stars 
(see Fig. |). 

Some specific observational complications arose, how- 
ever, since they are moving objects with respect to the 
sky background and show periodic variations of intensity 
due to rotation and varying distance from both the Earth 
and Sun. Selection criteria were: well understood rota- 
tional behaviour, small lightcurve amplitude, some form 
of independent size determination (either direct imaging 
or via occultation measurements), good visibility during 
the ISO mission, and availability of sufficient observational 
data from visible to submillimetre wavelengths. The se- 
lection resulted from a combination of extensive ground- 
based observing campaigns at thermal wavelengths from 
the IRTF, UKIRT, and JCMT, additional visible wave- 
length lightcurve measurements, and a series of FIR ob- 
servations from the now retired KAO. These observations 
were used to confirm the validity of the subsequent ther- 
mophysical modeling effort. 

Initially, modified versions of the Standard Thermal 
Model (Lebofsky|l98|) were used. For the final ISOPHOT 
photometric calibration, a thermophysical model (TPM) 
assuming a rotating ellipsoid and parameterising heat con- 
duction, surface roughness, and scattering in the regolith 
was adopted for 10 asteroids (Miiller & Lagerros 1998| ). 
The TPM is capable of producing thermal lightcurves 
and spectral energy distributions for any time, taking 
into account the real observing and illumination geom- 
etry. The overall comparison between model predictions 
and our large sample of MIR, FIR, and submillimetre ob- 
servations (about 700 individual measurements between 2 
and 2000 /im) demonstrated an accuracy of approximately 
±10% across a wide wavelength range from 10-500 /xm. 
In a few cases, existing direct measurements of shape sig- 
nificantly improved the model accuracies. 



In a recent study of the accuracy of the TPM, 
ISOPHOT observations of asteroids were independently 



calibrated in the FIR (Miiller & Lagerros |200l| , gOOlaj). In 
this work, parts of the FCS power curves were established 
using measurements only of planets or stars. These were 
used then to calibrate the asteroid observations and com- 
pare the results to the TPM predictions. It was shown 
that observations and predictions for 1 Ceres, 2 Pallas 
and 4 Vesta agree within 5 %. For 3 Juno, 10 Hygiea, 54 
Alexandra, 65 Cybele, and 532 Herculina the agreement is 
still within 10 to 15 %. For the 2 objects without indepen- 
dent ISO data the comparison to ground based thermal 
observations gave r.m.s. values of 14% for 313 Chaldaea 
and 29 % for 106 Dione. Note that the r.m.s. values in- 
clude also observational errors, which are dominated by 
the structured FIR sky background and calibration uncer- 
tainties. Modeling limitations are mainly due to uncertain- 
ties about the exact asteroidal shapes. The final modeling 
results were also confirmed to an accuracy of better than 
10 % via comparison with rotationally resolved target ob- 
servations from the IRAS database. Table |l| lists the solar 
system objects that were used for ISOPHOT calibration. 

5.4. In-Band Power Calculation 

To determine the responsivity of a detector we need to 
calculate the in-band power. Given the power distribution 
of a celestial standard in the focal plane, i.e. the normal- 
ized point spread function^, PSF{x,y, X), times the flux 
density of the source, F\(X) = {c/}?)F^, we define the in- 
band power as the integral over the product of i) the flux 
density of the source, ii) the PSF, iii) all spectral transmis- 
sions Tf{X) and iv) the relative spectral response function 
of the detector i?p(A), normalized to unity at its peak re- 
sponse. The integration is performed over all wavelengths 
and over the area of the detector aperture. 



Tf{X)Rp{X)FxiX) jJPSF{x, y, X)dxdy 



dA.(l) 



Wavelength-independent losses due to mirrors are ac- 
counted for by the factor T^. The indices p and / indicate 
dependence on pixel or filter, respectively. 

The filter transmissions (including the out-of-band re- 
jections over a wide wavelength range) and the relative 
response functions of the detectors were measured in the 
laboratory under "cold" conditions (at their anticipated 
operating temperatures in-orbit) with a Fourier-transform 
spectrometer (Schubert 1993). To calculate T^, we as- 
sumed the reflectance of each mirror surface to be 98 % 
and wavelength-independent. Thus we obtain = 0.98", 
with n being the number of mirrors in the optical paths, 
which is 6, 7, 7, 6 and 5 for the detectors PI, P2, P3, ClOO 
and C200, respectively. For PI an additional multiplica- 
tive factor of 0.95 was included to account for losses due 
to the Fabry lens. 



^ The integral of the PSF over the focal plane is normalized 
to 1. 
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We note that the in-band power is only a fraction of 
the total incident IR power, since the formula includes 
the wavelength-dependent relative response function of 
the detector material. Normalizing this response function 
to 1.0 at its peak allows us to define the detector respon- 
sivity as characteristic of the detector pixel, independent 
of the filter used. Thus Psrc represents only the fraction of 
the in-band power that actually initiates a photocurrcnt 
in the detector material. 

Since the wavelength dependence of the PSF within a 
filter band is small, we simplified Eq. to 



-Psrc — ^ /] 



PSFf 



/ [r/(A)i?p(A)^^;,(A)]dA. 



(2) 



by replacing the integral over the aperture and the PSF 
by two factors: (i) The area of the telescope aperture. A, 
and (ii) the PSF-correction factor, /psF/ „ , dependent on 
filter, /, and aperture, a. This second factor describes - for 
a given combination of filter and aperture - the fraction 
of the PSF that actually enters the detector aperture. 

For the initial calculations we used an analytical 
telescope model, taking into account the primary and 
secondary mirrors. Subsequently we replaced this by a 
more accurate model (Okumura 200C ) based on numer- 
ical Fourier transformation, that also includes effects of 
the support structure of the secondary mirror. The PSF 
factors were calculated as the ratio of the in-band power 
resulting from Eq. 0, integrated over the detector aperture 
and the total in-band power in the PSF, i.e. Eq. |l| inte- 
grated over an infinite aperture. A h'Fi, — const, spectrum 
was adopted as source SED. 

Thus far, our intent has been to determine the in-band 
fluxes of known standards to characterise the detectors 
and to calibrate the FCS. Once the FCS is calibrated, 
however, we can use the instrument to determine the in- 
band powers of other celestial sources. To convert these 
back to flux densities, we must invert Eq. ^. For regular 
target observations the source SED is generally unknown 
so we need a common reference spectrum with a defined 
colour. We adopt a spectrum of the form vF^, = const, as 
Beichman et al. ( |1988D did for IRAS. This enables us to 
replace Eq. by the much simpler form 



Psrc — F„ Clf /pSF/,a, 



(3) 



where CI / are tabulated factors for the "reference wave- 
length" (see Table ||) of each filter. The photometric re- 
sults, F^, provided by the off-line processing software 
(OLP) or the interactive analysis (PIA) to the observer are 
expressed in terms of flux densities at the proper reference 
wavelength of the filter used. The choice of the reference 
wavelength is arbitrary, as long as the reference spectrum 
is defined. For sources with colours deviating from the 
reference spectrum, a colour correction must be applied, 
that can be derived from Eq. |^ and |^. Tabulated colour 
correction factors were derived for a few typical spectra 
like blackbodies or power-laws. Table || shows the refer- 
ence wavelengths that are used to calculate ISOPHOT 
flux densities. 
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Table 2. The reference wavelengths, conversion factors, 
standard apertures and the calibrated flux-range for the 25 
ISOPHOT filter bands. The calibrated flux-range is where 
FCS 1 was directly calibrated against external standards. 
For revolutions prior to 94 sometimes smaller ranges ap- 
ply. Extrapolations are possible but may have larger abso- 
lute errors. For the C-subsection values are given for the 
central pixel of C'lOO and pixel 1 of C200. 



5.5. Observations 

The ISO observations to calibrate the internal sources 
were scheduled in blocks for each fllter band. Each block 
comprised measurements of a celestial source, celestial 
background, heated FCS and non-heated FCS. The mea- 
surements were obtained in different sequences, but in gen- 
eral the low-flux measurements were executed before the 
high-flux measurements to reduce the time constants of 
transients. These are shorter for transitions from low to 
high flux. The observations for the long- wavelength detec- 
tors P3, ClOO and C200 also included direct calibrations 
of the redundant FCS 2. For the other filter bands, FCS 2 
was calibrated relatively to FCS 1. 

To check the relative effective transmission of filter 
bands, measurement blocks of the same detector were fur- 
ther grouped together, so that responsivity drifts would 
be minimal during that time. We verified this by repeat- 
ing the first measurement of the sequence at the end. The 
results lead to correct ions to the filter transmissions, that 
are discussed in Sect. |6.2l 
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All detector pixels were calibrated individually. 
Therefore, measurements of the standard point sources 
were positioned at the centre of every pixel. To minimize 
the measurement time for the camera arrays, the space- 
craft performed a raster mode pointing sequence with a 
step size equal to the distance between pixel centres (46" 
for ClOO, 92" for C200). Larger raster dimensions (3x5 
for ClOO and 2 x 7 for C200) were chosen to obtain the 
background zero level also. Where no raster measurement 
was required, background measurements were performed 
in staring mode on positions within a 5' radius of the 
source. For aperture-photometry of faint sources, a one- 
dimensional raster of 3 points with the source at the centre 
position was scanned back and forth (nodding- mode), to 
enable elimination of the baseline drift (see Fig. ^) . 

The use of a raster observation on a fixed sky po- 
sition was also necessary because ISO could neither ob- 
serve positions with a relative offset from a tracked Solar 
System Object (SSO), nor perform raster observations 
while tracking. Therefore, all raster observations on SSOs 
had to be scheduled at fixed times with the telescope 
pointing adjusted to encounter the moving object. We 
made sure that the proper motion of the SSO due to 
the relative motion w.r.t. ISO was always below 60"/hr 
during observations. The contribution from ISO's motion 
relative to the centre of the Earth was only of the order 
of ±3"/hr, as science observations were made during the 
"slow" part of the orbit. During the maximum duration of 
a raster measurement on an SSO with ClOO (688 s), aster- 
oids would travel only « 12". Taking into account ISO's 
absolute pointing accuracy for fixed targets of about ±2", 
the worst case pointing error amounts to 6" -I- 2" = 8", 
which is still much smaller than the large detector aper- 
tures or pixel sizes used at long wavelengths. 

6. Photometric Corrections 

All corrections so far applied to the signal did not require 
quantitative knowledge of the IR flux illuminating the de- 
tector. However, in this section we will describe corrections 
that became necessary after a detailed photometric analy- 
sis of observations of known celestial calibration standards 
of various brightnesses. 

6.1. Signal Linearisation 

We defined earlier the responsivity, R, as a property of 
the detector that is expected to vary in time due to cosmic 
radiation and flux history, but not to depend on flux level. 
However, the actual behaviour of ISOPHOT detectors was 
found to be different. 

In its simplest form we calculate the responsivity from 
the signals measured on-source (S'src + S'bck) and off-source 
(S'bck) and the PSF-corrected in-band power, Pgrc, accord- 
ing to 



R ^ 



Q n. 
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Fig. 9. Naive plots of responsivities against total in-band 
power, assuming linearity (Eq. ^ for representative filter 
bands of all detectors. The total in-hand power includes the 
contributions from both source and background. Upper and 
lower diagrams show the responsivities before and after lin- 
earisation. The upper diagrams show a clear dependence 
on in-band power except for PS. The remaining scatter in 
the linearised data is due to the variation of the respon- 
sivity with time. 



where Cint is the capacity of the integrating amplifier. 
Plotting responsivities from measurements of various cal- 
ibration standards versus in-band power, however, shows 
strong correlations (see Fig. ^). Moreover these correla- 
tions are different even for filters of the same detector 
subsystem. 

If the measured responsivity depends in some way on 
the IR flux that falls onto the detector, R can also be 
expressed as a function of the photocurrent or the detector 
signal. We then rewrite Eq. ^ as 



R{Ssrc + S'bck) — 



(Ssrc + S'bck) Cii 



P 



bck 



(5) 



describing the responsivity of the detector when it is 
pointed at the position of the source, and further as 



R{Ss, 



Si 



bck) 



(Ssrc + Sbck) Cint 



(6) 



The index bck indicates the signal and in-band power 
at the background position. In our case, i?(S'bck) and 
R{Ssrc + S'bck) are not necessarily equal and the ratio of 
the signals measured on- and off-source is different from 
the ratio of the respective in-band powers. To determine 
R{S) requires knowledge about -R(Sbck) because of Eq. |[ 
Considerable errors are introduced in comparisons of very 
different signals if R is assumed to be independent of sig- 
nal. This applies in particular to multi-filter AOTs, multi- 
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Fig. 10. Derivation of filter-to-filter corrections. We display the data for some representative detectors. Each diagram 
shows the normalized responsivities for the indicated detector, derived from various filter sequences. Each sequence was 
measured on a calibration standard within a time interval, sufficiently short to avoid long term drifts. Each sequence 
is connected by a line. The standard deviation of the scatter of the relative filter measurements is 2 to 9 % for the 
P- detectors and 3 to 8 % for the C- arrays. 



aperture AOTs on extended sources, maps with large dy- 
namic ranges, and all FCS calibrations where the signals 
from FCS and sky are very different. For measurements on 
the sky, the biggest systematic errors are expected when 
the individual contributions of background and source 
alone are about equal. 

All further data reduction steps assume a linear sys- 
tem, so we decided to linearise the detector signals, S*, 
according to S' = T{S), where T is a continuously rising 
transfer function. The condition for T{S) is such that the 
results of S' inserted into Eq. |^ should no longer corre- 
late with the detector signal, but should be distributed 
around a constant responsivity with minimum dispersion. 
The absolute value of this constant is arbitrary because 
it cancels out in the flux calibration. We normalized it to 
the median of the nonlinearised responsivities calculated 
according to Eq. ^. It should be noted that application of 
the signal linearisation requires the prior subtraction of 
the dark signal. 

Although the FCS calibration data comprise a large 
number of measurements, the number of independent re- 
sponsivities distributed over the entire flux range per filter 
band ranges only between 14 and 28. We developed a man- 
ual fitting method to determine the linearisation tables, by 
displaying the results of Eq. ^ and interactively modifying 
the linearisation function, T{S). This approach led to a 
better exploration of parameter space, a smoother transi- 
tion at the boundaries, and a more tolerant treatment of 
outliers without the need to develop sophisticated fitting 
algorithms. We increased the objectivity of the method 
by keeping the shape of T simple, generally close to a 



polynomial of second or third order. We also minimised 
differences between pixels of the array detectors by fitting 
all pixels of one detector in parallel on the same screen. In 
spite of the large scatter of responsivities, simple functions 
T were found that significantly increase the consistency of 
data measured at different flux levels. 



Signal linearisation is performed by interpolation in 
a lookup table spanning the full range of possible sig- 
nals. This table is calibrated only in the range covered by 
calibration standards but, in practice, signals occur out- 
side this range that are higher, lower, or even negative in 
the case of noisy measurements close to the dark current. 
Reasonable entries for these cases had to be defined. For 
high signals beyond the point where data are available, 
we maintained a constant responsivity. Small signals are 
constrained so that an already dark-subtracted zero sig- 
nal is not changed by signal linearisation. Therefore, we 
continue T{S) between the smallest point in the lookup 
table that is still determined by valid data, and zero, by a 
straight line. This assumes that the responsivity does not 
vary in that range, c) The most difficult cases are nega- 
tive signals, because negative fluxes are simply undefined 
and appear only as a result of noise. The number of such 
cases actually encountered is low due to the reduced noise 
after signal averaging over periods of constant flux (SCP 
level). We decided to continue the transfer function such 
that T(5) = ~T{—S), to cover all situations that occur. 
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6.2. Filter-to-Filter Calibration 

Calculating the detector responsivity from standard 
source observations according to 



(Ss- 



'pj 



Shckpj ) C'int 



7 



(7) I 



after application of all signal corrections described so far, 
still yields different results for the same detector pixel, p, 
measured in different filters, /, during the same batch of 
observations. The ratios of responsivities between filters 
are reproducible for different calibration sources and can- 
not be explained by responsivity variations with time. 

The worst case is detector PI, where R varies by a 
factor of 4 between filters P_3.29 and P_16 (see Fig. p^ . 
Currently, the wavelength-dependent flux losses of this de- 
tector remain unexplained. We speculate, however, that 
the sapphire substrate that mechanically supports the de- 
tector crystal within the kidney-shaped aluminium cavity 
might play a role. Sapphire is transparent in the NIR and 
starts to absorb beyond about 6 /im onwards, in the mid- 



dle of Pf 's wavelength range (Tropf & Thomas |1998| ). The 
absorption peak is at wlT.S /im. Some fundamental lat- 
tice vibration modes occur in the reststrahlen region above 
13 /im. Top and bottom of both detector and substrate are 
gold-coated. The remaining effective area of the detector, 
where radiation from within the cavity can enter, is only 
twice that of the substrate. Given the high refractive index 
of Si, around 3.4, and the consequently higher reflectivity 
compared to AI2O3 (refractive index of only 1.7), we esti- 
mate that about the same number of photons enters both 
detector and substrate. The cavity was designed to coun- 
teract the high refractive index of the detector material by 
maximizing the number of reflections within it. Assuming 
that the sapphire substrate turns "dark" at longer wave- 
lengths, it is conceivable that, with every reflection, a sub- 
stantial fraction of the radiation enters the substrate and 
is attenuated. The effective response of the detector sys- 
tem would thus deteriorate. Another IR absorber present 
within the cavity is glue (Stycast). Unfortunately the ques- 
tion of Pi's quantum efficiency remains largely academic, 
so that the effort of a much more detailed analysis is not 
justified. 

The other detectors show less dramatic variations by 
factors up to 2. Moreover, the ratios found between differ- 
ent filters are not the same for different pixels of the same 
C-detector array. We attribute this to a projection of spa- 
tial nonuniformities of the filter surface onto the detector 
array, because the ISOPHOT filters are not located at the 
pupil of the optical path, but close to the detectors. The 
overall variation between filters of all detectors, except PI, 
is most likely due to optical misalignments and diffraction 
effects that are neglected in our ideal instrument model. 

We applied a correction to the responsivity by intro- 
ducing a matrix, XpJ- If works equally for P- and C- 
detectors. To obtain a wavelength independent detector 
responsivity, we split R' of Eq. |^ into filter-independent 
and filter-dependent parts, so that RpXpj = ^pf- "^^^ 




1.0 10.0 
FCS heating power [mW] 



Fig. 11. Constructing a power curve from calibration 
data. The names indicate the sources that have been used 
to calibrate the respective data points. 



matrix, Xpji was established from multi- filter measure- 
ments on calibration standards, executed so that the 
change of responsivity with time could be neglected for 
a filter sequence. After eliminating poor quality measure- 
ments, the responsivities measured for a detector pixel, 
p, in a filter sequence were renormalized. The factors per 
sequence were chosen to achieve the best match among 
all sequences for a given pixel. After calculating the av- 
erage responsivities, Rpj, per filter, /, we normalized 
to the maximum responsivity found for that pixel using 
XpJ = -Rp,//inax(i?pj), thus assuming that the matrix 
describes flux losses w.r.t. the ideal model. The normalized 
data and the resulting filter-to-filter corrections, Xp,/, for 
some detector pixels and filter bands, are shown in Fig. ^ 
The standard deviation of the scatter of the relative filter 
measurements is 3 to 8 % for the C-arrays and 2 to 9 % 
for the P-detectors. 

7. Derivation of FCS power curves 

The FCSs were calibrated against celestial standards in 
all 25 filter bands. These measurements yield the relation- 
ship between in-band power, P, and FCS heating power, 
/i, for all combinations of filter band, /, detector pixel, p, 
and FCS. Signals are measured differentially. The source 
signal is obtained from the difference between on-source, 
S'src, and off-source, S'bck, signals while the FCS signal 
is calculated as the difference between the signals mea- 
sured with heated FCS(S'fcs) and cold FCS(S'str)- Since 
the measurements of FCS and celestial standard were per- 
formed in the same filter, the responsivity does not need 
to be explicitly calculated. We used 



p/ /, s, _ ('S'fCs(^) - 'S'str) ^src 



>5's] 



Sh 



ck 



(8) 



where Ppcgih) is the in-band power incident on the de- 
tector from the FCS and Pgrc is defined in Eq. |2[ Fig. |l| 
shows an example of the power curve resulting from ob- 
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♦ PIA pixel numbering 
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Fig. 12. Relative illumination of the array detectors ClOO 
and C200 (illumination matrices) by FCS 1 (valid after 
rev. 93). Pixel numbers and ISO coordinate axes are given. 



treatment of these cases is still poorly understood and is 
beyond the scope of this paper. For use with standard and 
very similar apertures, the in-band powers of the PHT-P 
power curves were divided by the area of the standard 
aperture of the measurement, and are expressed in units 
of W/mm^. 



7.2. Extrapolation 

During the mission, FCS heating powers were used that 
fell beyond the calibrated range of heating powers. We 
used the grey body model described in Sect. || to extrap- 
olate the power curves. Including the spectral transmis- 
sion of the beam splitter of the FCS and replacing the 
celestial source SED by the grey body spectrum, the in- 
band powers are calculated in accordance with Eq. ^ The 
free parameters a and (3 and the effective solid angles of 
the detector subsystems are obtained from least-squares 
fits to the FCS curves (Schulz 1993 ). Comparison of the 
model to the in-orbit data shows good consistency, espe- 
cially at longer wavelengths. At shorter wavelengths the 
consistency is poorer. Extending the measured FCS power 
curves by the models provided the values outside the range 
of measured heating powers. This range is flagged in the 
associated calibration tables. 



servations of different calibrators for pixel 4 in the C_105 
filter band. 



7.1. Illumination Matrices and FCS Apertures 

FCS power curves for different pixels of the same C- 
array differ by constant factors. These factors result from 
both the nonuniform illumination of the detector arrays 
by the FCS and the putative spatial inhomogeneities in 
the filter. For each filter, these factors constitute the il- 
lumination matrix, Ppj, as shown in Fig. |l^. Thus the 
FCS curve for each filter is split into a pixel-independent 
part, PFCSfih)i and a pixel-dependent factor, F, so that 
Ppcs fW — ^pJ PpcSfih)- The FCS curve is normahzed 
such that Fpj, averaged over all pixels becomes 1, for a 
given filter. The resulting FCS power curves for FCS 1 are 
shown in Figs |l^ and |lj. Between revolution 93 and 94 
the TRS 2 of FCS 1 , which was routinely used with the 
long-wavelength detectors, abruptly increased its bright- 
ness by about a factor of 2. This change necessitated a re- 
calibration of the corresponding filter bands and resulted 
in a second set of FCS calibration tables for P3, ClOO 
and C200. Fig. [| shows the FCS power curves that ap- 
ply to the period before revolution 94. The non-uniformity 
of the FCS illumination also affects the FCS signal mea- 
sured in different apertures of the P-detectors. For a fixed 
heating power the FCS signal is not proportional to the 
area of the aperture. Strong deviations exist for apertures 
much larger or smaller than the standard apertures. The 



73. FCS Straylight 

Ground calibrations already showed an increased signal 
compared to the dark signal when measuring the signal 
from a "cold" FCS. Considering the internal ISOPHOT 
geometry, this additional signal is most likely picked up 
from ra diatio n within the field of view of the telescope 
(Schulz |l993 ). The in-orbit observations, where the cold 
FCS signal was found to be correlated with the back- 
ground signal, confirm the presence of this kind of stray- 
light for C200 (CornwaU |l997| ). ClOO and P3 show no de- 
tectable straylight component. P2 shows a constant value 
of 0.02-0.04 V/s for a signal > 0.2 V/s (on source) and a 
steep fall-off for lower signals. Straylight in PI is always 
below 0.07 V/s but shows a nonlinear dependence on sig- 
nal. The nonlinear dependences for both detectors PI and 
P2 are quite uncertain but appear to be due to the varying 
dark signal and systematic errors in the transient fit at low 
signals, rather than to straylight, which is expected to be 
proportional to the intensity of the incident flux compo- 
nents. We conclude, therefore, that FCS straylight should 
be taken into account for C200, but can be neglected for 
the other cases. We determined the straylight at the FCS 
position measured with C200 to be 5 % of the background 
level. 
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Fig. 13. FCS power curves of all filter bands. For PS, ClOO and C200 we show the versions valid after the change 
of FCS 1/TRS 2 in rev. 94. The interpolated parts are indicated by solid or dashed lines with small squares at the 
boundaries, while the model extrapolations are indicated by dotted lines. 




0.1 1.0 10.0 0.1 1.0 10.0 0.1 1.0 10.0 

Healing Pmver |mW| Heatini: Power |m\V| Healing Power |mW| 



Fig. 14. FCS power curves of all filter bands of the long wavelength detectors before the change in FCS 1/TRS 2. 



8. Flux Calibration 

We can now detail the procedure to calibrate ISOPHOT 
observations. First, the detector responsivity is deter- 
mined from the FCS measurement using 



(9) 



^^08/(^)1 Tp,/ and Xpj are taken from the calibration ta- 
bles. The FCS straylight signal Sstij, j was actually mea- 
sured only in absolute photometry AOTs and in dedicated 
FCS calibration measurements. In all other cases it is re- 
placed by the time-dependent dark signal (see Sects. 4.7 
and 7.3). C200 should be corrected for the 5% stray- 
light contribution of the background. Second, the in-band 



B. Schulz et al.: ISOPHOT - Photometric Calibration of Point Sources 



17 



power, Pf, for all filters of the detector subsystem is cal- 
culated as 



P, 



SsTCp^f C'int 



(10) 



The transformation to a flux density, i^^, at the reference 
wavelength of the filter (e.g. 12 /im for the P_11.5 filter) 
is done by changing Eq. to 



CI/ /PSF/,, 



(11) 



The flux density of the point source is obtained by sub- 
tracting the flux density derived from a measurement on 
the background, following the same calibration steps as 
for the source observation. The resulting flux density is 
calculated at the reference wavelength of the filter for an 
assumed vF^ — const, spectrum and still has to be colour- 
corrected in accord with the actual source spectrum. It 
should be noted that, when the point source is not placed 
at the centre of the detector pixel or aperture, a different 
PSF correction factor must be applied. In particular, in 
observing modes with C-detectors, the sum of the fluxes 
found in all pixels can be calculated and a PSF correction 
factor, fpsFf „ , determined for the entire array, is applied. 

8.1. Accuracy 

We assess the overall consistency of the point source cali- 
bration in standard apertures by comparing the measured 
FCS in-band fluxes for all available heating powers with 
the FCS power curves. The scatter is representative of 
the error budget, except for any systematic biases in the 
SEDs of the calibration sources. Each measurement passed 
through the standard signal processing chain and thus 
contains the same uncertainties as the science observa- 
tions. The additional automatic transient fit we applied 
during the calibration data processing is equivalent to im- 
provements the general observer can achieve through man- 
ual processing by Interactive Analysis. We plot the ratio 
of the in-band power taken from the calibration tables for 
the FCS (FCS power curves times FCS illumination ta- 
bles) and the measured FCS in-band powers. To indicate 
the flux regimes, we plot the ratios versus in-band power 
on the detector. All data from a single subinstrument ap- 
pear together in one diagram. Figs |l^ and ^ show the 
ratios derived for the bulk of the mission for FCS 1 and 
the ratios that apply for FCS1/TRS2 before revolution 
94. As expected, the consistency improves with flux (i.e. 
the scatter diminishes), because statistical noise becomes 
less important and time constants of transients decrease. 

To obtain another quantitative representation, we con- 
sidered the deviations for individual fllters and computed 
the maxima and the r.m.s. We also calculated the r.m.s. 
of the portions left and right of the dotted lines in Figs ^ 
and |l6|, to seek any differences between higher and lower 
fluxes. Table || contains the values applicable for most of 
the mission while Table ^ shows the values for the time 
before rev. 94. These flgures and tables give an overview of 
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Table 3. Comparison of derived FCS power curves with 
measured calibration data. The deviations are given in per- 
cent, n is the number of measurements used, max is the 
maximum deviation found and a is the r.m.s. deviation. 
The last two columns a\o and dhi give the r.m.s. deviations 
in the low and high flux regions separated by the vertical 
dashed lines in Fig. 15 and\ld^ 



the accuracy expected for an observation of a point source 
in the standard aperture, including measurements of back- 
ground, source, FCS and cold FCS. Since the FCS stray- 
light measurement can in most cases be replaced by the 
dark signal, the validity of this assessment extends to even 
a larger number of observations. For multi- filter measure- 
ments, where an FCS measurement in only a single filter 
is available to determine R (Eq. this FCS calibration 
can be transferred to a sky measurement in another filter 
of the same subsystem, because the ratios of transmis- 
sions between filters have been determined (see Sect. 3.2). 
However, this transfer also accretes the uncertainty of the 
factor Xp.f (~ 10%) for any filter different from the one 
that was actually used during the FCS measurement. 

For a single staring observation performed in a stan- 
dard aperture, we find scatters in the range from 3 to 29 %, 
depending on detector, filter band and flux. The inhomo- 
geneous distribution of data points suggests that system- 
atic effects still dominate the error budget. Since the data 
points are few and the properties of the distribution, which 
may not be Gaussian, also depend on flux or signal level, 
the results must be interpreted with care. The average 
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Fig. 15. Calibration accuracy derived from calibration consistency. The diagrams show the ratio of the FCS power 
curves times the illumination matrices and the data points derived for the P-detectors, plotted against FCS in-band 
power. The residuals of all filters of a subsystem are combined within one plot. The scatter is a measure of the error 
budget for single observations. The vertical dashed line separates the low and high flux intervals. 
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22 


25 


19 


C_70 


108 


29 


14 


13 


14 


C_90 


109 


31 


17 


20 


14 


C_100 


94 


39 


16 


18 


12 


C.105 


92 


41 


17 


16 


17 


C.120 


49 


17 


7 


7 


6 


C.135 


41 


16 


6 


6 


7 


C_160 


58 


23 


8 


8 


8 


C.180 


35 


13 


5 


4 


6 


C_200 


27 


15 


7 


7 





Table 4. Same as previous table but for the period before 



scatter (10%) over all filters of detector PI appears to be 
the best metric for the accuracy that can be expected, but 



it is likely to be better. A similar result applies to detector 
P2. Only P3 shows a much larger scatter of about 20% 
as a consequence of its degraded performance in the ra- 
diation environment in space compared to the laboratory 



(Lemke et al. 1996). Better statistics are available for the 
9 pixels of ClOO, that perform better than P3 although 
the detector material (Ge:Ga) is the same. We find an av- 
erage scatter of 13%. A better background subtraction, 
due to the raster mapping method we applied for the C- 
detectors, may also contribute to the reduced scatter. The 
results for the C200 detector also provide good statistics. 
The average scatter is about 10%, fairly consistent over 
all filters. 

These numbers serve as a guideline for the astronomer 
as to what accuracies to expect for a staring or raster 
observation that underwent a complete data reduction, 
involving all reduction steps and all necessary support- 
ing measurements, i.e. background, FCS, FCS straylight. 
To assess the uncertainty of the FCS power curves alone, 
however, they only represent upper limits, since statisti- 
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Fig. 16. Calibration accuracy derived from calibration consistency. The diagrams show the ratio of the FCS power 
curves and the data points derived for the C-detectors, plotted against FCS in-band power. The residuals of all filters 
of a subsystem are combined within one plot. The scatter is a measure of the error budget for single observations. The 
vertical dashed line indicates the separation between low and high flux intervals. 



cal uncertainties decrease with the number of data points 
used. 

The uncertainties obtained are similar to or larger than 
the accuracies quoted for the SEDs of the celestial stan- 
dards. This limits our ability to discriminate problematic 
standards in comparison to others. A first analysis is given 
by Schulz (|200lD. 



8.2. Long-term Monitoring 

In order to test the stability and the reproducibility of 
the system, a set of non-variable celestial sources was 
monitored at regular intervals throughout the mission. A 
bright and a faint source were repeatedly measured in the 
standard filter band of each detector, using the standard 
apertures in the case of the P-detectors. For the bright 
source we used the planetary nebula NGC6543, which 
was the prime photometric reference for IRAS due to its 
good visibility a nd de tectability from the MIR to the FIR 
(Beichman et al. 1988 ). The faint sources were chosen from 



the list of stellar calibration standards with good visibility 
throughout the mission (see Table I). 

The measurements on NGC6543 were repeated every 
5 weeks on average; the fainter sources were measured ev- 
ery 2 weeks. The planetary nebula was measured at the 
centre of each detector system, after performing a back- 
ground measurement at a position about 5' away. The 
faint sources were observed either in nodding mode for 
the P-detectors or with a small raster of 3 x 3 or 4 x 2 for 
ClOO and C200 respectively, to minimize the uncertainties 
due to background subtraction. We find the reproducibil- 
ity of measurements on bright sources to be around 2- 
3%, except for PI and P3 for which this is 6-7%. The 
fainter sources show larger scatter because uncertainties 
of source, background, FCS and FCS straylight contribute 
equally to the error, compared to the uncertainties only of 
source and FCS that contribute at higher fluxes. Photon 
noise and glitch noise also contribute relatively more at 
low fluxes. The reproducibility below 1 Jy ranges from 7 
to 12%. A summary is given in Table |[ 
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source 


filter 


apert. 






a 








[Jy] 


[Jy] 


[%] 


HD172323 


P_11.5 


52" 


0.085 


0.087 


10.0 


HR5986 


P_25 


79" 


0.667 


0.588 


12.3 


HR7310 


P.60 


180" 


0.667 


0.574 


7.0 


HR7310 


C.lOO 


- 


0.229 


0.250 


10.6 


HR6705 


C.160 




0.672 


0.631 


7.2 


NGC6543 


P_11.5 


52" 




6.8 


5.9 


NGC6543 


P_25 


79" 




107 


1.8 


NGC6543 


P.60 


180" 




144 


6.6 


NGC6543 


C.lOO 






59.5 


2.0 


NGC6543 


C.160 






15.9 


2.6 



Table 5. Overview of results of the monitoring programme 
of faint and bright flux standards. The sources, filters 
and apertures used appear in the first three columns; the 
fourth and fifth show the flux densities calculated from 
model SEDs and the observed error-weighted mean values. 
Flux densities are not colour- corrected, i.e. they assume a 
vF^ ~ const, spectrum. The last column shows the stan- 
dard deviation of the distribution of measurements. 



9. Summary 

We have described the photometric calibration of 
ISOPHOT using celestial standards. We used point 
sources either stars, planets or asteroids, to cover the full 
wavelength range. Standards at all flux levels were used 
to calibrate nonlinearities. 

Because the detector responsivity varied with time, 
stable internal IR sources were calibrated against celes- 
tial standards and measurements of these references were 
included in each scientific observation. The reproducibil- 
ity of results derived from the ratios of these signals was 
verified for high and low flux levels. We corrected for non- 
ideal effects like nonlinearities, imperfect optics and de- 
tector transients, and reduced their influence on the mea- 
surements to a minimum. 

The photometric calibration of the internal sources ini- 
tiated the definition of empirical FCS power curves and 
FCS illumination tables as well as the introduction of 
matrices of calibration factors, to account for pixel- and 
filter-dependent attenuations, which are probably caused 
by misalignments, diffraction and spatial gradients in fil- 
ter transmissions. A major step in achieving a consistent 
picture was the correction for the flux dependence of the 
responsivity. A set of formulae describing the path from 
raw detector signals to flux densities was derived, which 
uses the detector responsivity to characterise the state of 
a detector pixel at the time of the observation. 

To quantify the accuracies achievable by absolute pho- 
tometry, a comparison was made of the derived calibra- 
tion curves of the internal sources and the measured data. 
Depending mainly on flux level, the achievable uncertain- 
ties are around 10 to 20%, but can exceed 30% under 
exceptional conditions and at low flux levels. The fllters 



of each subsystem were calibrated relative to each other 
with accuracies better than 10%. 

The use of point sources and standard apertures for 
each filter for the aperture photometer defines a base- 
line within the parameter space, which serves as a refer- 
ence for the calibration of further instrument modes like 
chopped measurements, multi-aperture measurements and 
mapping. 
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